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1 Introduction 

The purpose of this paper is to describe two algorithms for computing with 
word-hyperbolic groups. Both of them have been implemented in the second 
author's package KBMAG ||Holt, 1995| . 



The first is a method of verifying that a group defined by a given fi- 
nite presentation is word-hyperbolic, using a criterion proved by Papasoglu 
in 



papasoglu, 1994|1 , which states that all geodesic triangles in the Cayley 
graph of a group are thin if and only if all geodesic bigons are thin. This is 
very similar to an algorithm described in [[Wakefield, 1997|| , but it contains a 



simplification which appears to improve the performance substantially. It also 
improves a less developed approach to the problem described in ||Holt, 1996 



The second algorithm provides a method of estimating the constant of 
hyperbolicity of the group with respect to the given generating set. We do 
not know of any previously proposed general method for solving this problem 
that has any prospect of being practical. Our current implementation is 
experimental, and is very heavy on its use of memory for all but the most 
straightforward examples, but it does at least succeed on examples like the 
Von-Dyck triangle groups and the two-dimensional surface groups. 

Both of them follow a general philosophy of group-theoretical algorithms 
that construct finite state automata. This approach was originally proposed 
in the algorithms described in Chapter 5 of [Epstein et ai., 1992|| for comput- 



ing automatic structures, and employed in their implementations for short-lex 



structures described in [[Holt, 1996| . The basic idea is first to find a method of 
constructing likely candidates for the required automata, which we shall call 
the working automata, The second step is to construct other (usually larger 
and more complicated) test automata of which the sole purpose is to verify 
the correctness of the working automata. In the case when this verification 
fails, it should be possible to use words in the language of the test automata 
to construct improved versions of the working automata. One practical dif- 
ficulty with this approach is that experience shows that incorrect working 
automata and the resulting test automata are much larger than the correct 
ones, so it can be extremely important to find good candidates on the first 
pass. 
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The two algorithms dealt with in this paper are described in Sections ||| 
and |j. Some details of their performance on some examples are presented 
in Section |5|. Unfortunately, we need to recall quite a lot of notation from 
related earlier works, and we do this in Section ^|. 

There are a number of computational problems in which it is useful to 
know hyperbolic constants which are either the same as those in this paper, 
or are closely related to them. Some examples of such problems follow. 

In another paper we will describe a linear algorithm to put a word in a 
word-hyperbolic group into short-lex normal form. The linear estimate is due 
to Mike Shapiro (unpublished). Our method, which is a bit different, may 
be usable in practice, though the ideas have not yet been implemented. The 
standard algorithm for converting a word into normal form in an automatic 
group is quadratic, as shown in [[Epstein et al., 1992] . 

We also plan to show how to construct the automata which accept 1) all 
bi-infmite geodesies, 2) all pairs of asymptotic geodesies and 3) all pairs of bi- 
infinite geodesies which are within a finite Hausdorff distance of each other. 
Such algorithms are necessary if one is to have any hope of a constructive 
description of the limit space of a word-hyperbolic group, starting with gen- 
erators and relations of the group. Some of these automata are also needed 
in Epstein's nlog(n) solution of the conjugacy problem (not yet published). 



2 Notation 

Throughout the paper, G will denote a group with a given finite generating 
set X. The identity element of G will be denoted by 1g- Let A = X U X" 1 , 
and let A* be the set of all words in A. For m,d6 A*, we denote the image 
of u in G by u, and u =g v will mean the same as u = v. For a word u G A* , 
l(u) will denote the length of u and u(i) will denote the prefix of u of length 
i, with u{i) = u for i > l{u). 

Let T = Tx{G) be the Cayley graph of G with respect to X. We make 
T into a metric space in the standard manner, by letting all edges have unit 
length, and defining the distance d(x, y) between any two points of T to be 
the minimum length of paths connecting them. (The points of T include 
both the vertices, and points on the edges of T.) This makes T into a geodesic 
space, which means that for any x, y G T there exist geodesies (i.e. shortest 
paths) between x and y. For g G G, 1(g) will denote the length of a geodesic 
path from the base vertex 1^ of T to g. 

A geodesic triangle in T consists of three not necessarily distinct points 
a,b,c together with three directed geodesic paths u,v,w joining be, ca and 
ab, respectively. The vertices a, b, c of the triangle are not necessarily vertices 
of T; they might lie in the interior of an edge of T. 

There are several equivalent definitions of word-hyperbolicity. The most 
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convenient for us is the following. Let A be a geodesic triangle in T with 
vertices a,b,c and sides u, v, w as above. (Hence l{u) = d(b,c), etc.) Let 
p(a) = (l(v) + l(w) — l{u))/2 and define p(6),p(c) correspondingly. Then 
p{b) + p(c) = l(u), so any point d on u satisfies either d(d,b) < p{b) or 
d(d,c) < p(c), and similarly for v and w. The points d,e,f on u,v,w with 
d(d, b) = p{b) and d(d, c) = p(c), etc., are known as the meeting points of the 
triangle. 

In a constant curvature geometry (the euclidean plane, the hyperbolic 
plane or the sphere), the meeting points of a triangle are the points where 
the inscribed circle meets the edges. In more general spaces, such as Cayley 
graphs, the term inscribed circle has no meaning, but the meeting points can 
still be defined. 




Figure 1. This picture shows the meeting points as the intersection of the in- 
scribed circle with the edges of the triangle in the case of constant curvature 
geometry. 

Suppose a, b and c are vertices in the Cayley graph. Then the meeting 
points are also vertices if and only if the perimeter l(u) + l(v) + l(w) is even. 

Let 5 G Then we say that A is 5 -thin if, for any r e R with < 
r < p(x), the points p and q on v and w with d(p,a) = d(q,a) = r satisfy 
d{p, q) < 5, and similarly for the points within distance p(b) of b and p(c) of 
c. We call such points p and q A- companions, or, if A is understood, just 
companions. Note that the definition makes sense even when the triangle 
A is not geodesic — we measure distances along the edges of the triangle. 
Normally companions are distinct, but there can be many situations where 
they coincide — for example two geodesies sides of a triangle could have an 
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intersection consisting of a disjoint union of three intervals. Mostly points 
on the triangle have exactly one companion, but the meeting points normally 
have two companions — once again, in degenerate situations two or all three 
of the meeting points may coincide. 

The group G is called word-hyperbolic if there exists a 5 such that all 
geodesic triangles in F are 5-thin. (It turns out that this definition is in- 
dependent of the generating set X of G, although the minimal value of 5 
does depend on X.) The multi-author article [[Alonso et al, 1991f is a good 
reference for the basic properties of word-hyperbolic groups. 

We also need to recall some terminology concerning finite state automata. 
This has been chosen to comply with that used in pipit, 1996] as far as 



possible. The reader should consult [ [Epstein et al., 1992|| for more details on 
the definitions and basic results relevant to the use of finite state automata 
in combinatorial group theory. 

Let W be a finite state automaton with input alphabet A. We denote 
the set of states of W by S(W) and the set of initial and accepting states by 
T{W) and ^(W) respectively. In a non-deterministic automaton there may be 
more than one transition with a given source and label, and some transitions, 
known as e-transitions, may have no label. In a deterministic automaton, 
there are no e-transitions, at most one transition with given source and label, 
and W has only one initial state. (This type of automaton is named partially 
deterministic in [[Epstein et al., 1992|| .) In this case, we denote the unique 
initial state by cr (X) and, for each x G A and o G S(X), we denote the 
target of a transition from a with label x by a x if it exists. We can then 
define a u , for a G S(X) and u G A* , in the obvious way whenever all of the 
required transitions exist. The automata in this paper can be assumed to be 
deterministic unless otherwise stated. 

The automata that we consider may be one-variable or two-variable. In 
the latter case, two words u and v in A* are read simultaneously, and at 
the same speed. This creates a technical problem if u and v do not have 
the same length. To get round this, we introduce an extra symbol $, which 
maps onto the identity element of G, and let A' = A U {$}. Then if (u, v) is 
an ordered pair of words in A, we adjoin a sequence of $'s to the end of the 
shorter of u and v if necessary, to make them both have the same length. The 
resulting pair will be denoted by {u,v)\ and can be regarded as an element 
of (A* x A^)* . Such a pair has the property that the symbol $ occurs in at 
most one of u and v, and only at the end of that word, and it is known as a 
padded pair. We shall assume from now on, without further comment, that 
all of the two-variable automata that arise have input language A* x A* and 
accept only padded pairs. 

Note that if M is a two-variable automaton, then we can form a non- 
deterministic automaton with language equal to {u | 3t> : (u, v)^ G L(M)} sim- 
ply by replacing the label (x, y) of each transition in M by x. (This results 
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in an e-transition in the $.) We call this technique quantifying over 

the second variable of M. 

Following ||Holt, 1996| , we call a two- variable automaton M a word- difference 
automaton for the group G, if there is a function a : S(M) — ► G such that 

(i) a(a (M)) = 1 G , and 

(ii) for all x,y G and r G S(M) such that r^'^ is defined, we have 
oc(j( x ' yS >) = x~ x a{r)y. 

We shall assume that all states r in a word- difference automaton M are 
accessible; that is, there exist words u, v in A* such that a (M)^ = r. It 
follows from properties (i) and (ii) that a(r) = u _1 t>, and so the map a is 
determined by the transitions of M. 

Conversely, given a subset T> of G containing 1q, we can construct a word- 
difference automaton D with S(D) = V, o~q(D) = 1q and, for d, e G V and 
x,y E A^ a. transition d — > e with label (x, y) whenever x~ l dy —g e. The map 
a is the identity map. We call this the word-difference automaton associated 
with V. (We have chosen not to specify the set A(D) of accepting states 
of D, because this may depend on the context.) Having constructed the 
automaton, we throw away the elements of T> which are not accessible from 
the initial state o~o(D). 

If u, v G A*, then we call the set V = {u(i) v(i) \ i G Z, i > 0} the set 
of word- differences arising from (u,v). Then (u,v) is in the language of the 
associated word-difference machine D provided that vT x v G A(D). 

The group G is said to be automatic (with respect to X), if it has an 
automatic structure. This consists of a collection of finite state automata. 
The first of these, denoted by W, is called the word-acceptor. It has input 
alphabet A, and accepts at least one word in A mapping onto each g G G. 
The remaining automata M x , are called the multipliers. There is one of these 
for each generator x G A, and also one for x = 1q- These are two-variable, 
and accept (wi,w 2 )^ for Wi,w 2 G A*, if and only if Wi,w 2 G L{W) and 
w\x =g w 2 . See [[Epstein et al., 1992|| for an exposition of the basic properties 
of automatic groups. It is proved in Theorem 2.3.4 of that book that there 
is a natural construction of the multipliers M x of an automatic structure as 
word- difference machines. 

Now we fix a total order on the alphabet A. The automatic structure 
is called short-lex if the language L(W) of the word-acceptor consists of the 
short-lex least representatives of each element g G G\ that is the lexicograph- 
ically least among the shortest words in A* that map onto g. The existence 
of such a structure for a given group G depends in general on the generating 
set X of G, but word-hyperbolic groups are known to be short-lex automatic 
for any choice of generators. (This is Theorem 3.4.5 of [[Epstein et al, 1992|| .) 

The group G is called strongly geodesically automatic with respect to X if 
there is an automatic structure in which L(W) is the set of all geodesic words 
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from 1q to g for g G G. It is proved in Corollary 2.3 of pPapasoglu, 1994(1 that 
this is the case if and only if G is word-hyperbolic with respect to X (from 
which it follows that this property is independent of X). This result will be 
the basis of our test for verification of word-hyperbolicity in Section [| the 
procedure that we describe will verify strong geodesic automaticity. 

We shall assume throughout the paper that our group G = (X) is short- 
lex automatic with respect to X, and that we have already computed the 
corresponding short-lex automatic structure {W, M x \x G A U {1g}}- We 
assume also that the set Dm of all word- differences that arise in the multipliers 
M x together with the associated word- difference machine Dm in which 1q is 
the unique accepting state has been computed. These automata can be used 
to reduce (in quadratic time) words u G A* to their short-lex equivalent word 
in G and so, in particular, we can solve the word problem efficiently in G. 
The above computations can all be carried out using the KBMAG package 
||Holt, 199511 . 



3 Verifying hyperbolicity 

Papasoglu (| )Papasoglu, 1994]] ) has shown that a necessary and sufficient con- 
dition for a group to be word-hyperbolic is as follows. Let V be the Cayley 
graph with respect to some set of generators. The condition is that there 
is a number cp, such that, for any two geodesic paths u, v : [0,£] — > T 
parametrised by arclength, if u(0) = v(0) and u(£) = v(£), then, for all t 
satisfying < t < £, dr(u(t),v(t)) < cp. The least possible value of cp 
is called Papasoglu's constant. Such a configuration of u and v is called a 
geodesic bigon. In order to know that all geodesic bigons have uniformly 
bounded width, there is no loss of generality in restricting to the case where 
u(0) = v (0) is a vertex of the Cayley graph. We can also restrict to the case 
where u(£) = v(£) is either a vertex of the Cayley graph or the midpoint of 
an edge. It is unknown whether the uniform thinness of such geodesic bigons 
follows from the uniform thinness of the more special geodesic bigons with 
both ends vertices. 

Our algorithm not only verifies word-hyperbolicity for a given group, but 
also gives a precise computation of the smallest possible value of Papasoglu's 
constant. In fact, it gives even more precise information, namely the set of 
word- differences u(i) v(i), where u and v vary over all geodesic bigons with 
u(0) = v(0) a vertex and % varies over all positive integers. In all the examples 
we have looked at, we have observed that the number of such group elements 
is very much smaller than the number of group elements of length at most 
cp. This is important in practical computations. 

The algorithm for verifying word-hyperbolicity proceeds by constructing 
sequences WD n , GE n , GW n , T n , for (n > 0), of finite state automata. (These 
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letters stand respectively for 'word-difference', 'geodesic-equality', 'geodesic 
word acceptor' and 'test'.) In general, for n > 0, WD n will be a set of ele- 
ments of G containing {1g}, and WD n will be the associated word-difference 
machine in which 1q is the only accepting state. 

We define WT>i to be the set V M defined at the end of Section §, and let 
W be the short-lex word-acceptor. Then, when n > 1, we define GE n , GW n 
and T n as follows. We define the language 

L(GE n ) = {(u,v) G A* x A* | (u,v) G L(WD n ),v G L(W),l(u) = l(v)}. 

Recall that all elements of L(W) are geodesies and that the only accept state 
of WD n is 1q- It follows that in the previous definition, u = v and that u 
and v are both geodesies. Now we define the language 

L(GW n ) = {u G A* | 3v G A* : (u, v) G L(GE n )}. 

Again, u must be a geodesic. 
Then we define the language 

L(T n ) = {weA*\ L(GW n ) | 

3u : (w, u) G L{WD n ),u G L(GW n ), l(u) = l(w)}. 

Again, u and w are both geodesies in the previous definition. 




Figure 2. This illustrates the geodesic paths u, v and w described in Section |[ 

If L(T n ) is empty for some n, then the procedure halts. Otherwise, we 
find a geodesic word w G L(T n ), reduce it to its short-lex least representative 
v, and define WV n+ i to be the union of WV n and the set of word-differences 
arising from (w, v). Then we can define the automaton WD n+ i and construct 
the other automata for the next value of n. 

Theorem 3.1 The above procedure halts if and only if G is strongly geodesi- 
cally automatic with respect to X . 

PROOF: First note that, if L(T n ) is non-empty for some n, and contains the 
word w reducing to v G L(W), then the word-differences arising from (w,v) 
cannot all lie in WV n . For otherwise we would have (w,v) G L(WD n ), and 
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hence (w,v) G L(GE n ) and w G L(GW n ). But w has been chosen so that 
it is not in L{GW n ). Hence WV n+ i strictly contains WV n . Thus, if the 
procedure does not halt, then the number of word-differences arising from 
pairs of geodesies (u, v) with u =g v cannot be finite, and so by Theorem 
2.3.5 of ||Epstein et al., 1992| G cannot be strongly geodesically automatic. 



Conversely, suppose that the procedure does halt, and that L(T n ) is empty 
for some n > 1. We claim that L{GW n ) is equal to the set of all geodesic 
words. If we can show this, then it would follow from the definition that 
L(GE n ) is equal to the set of pairs of geodesic words {u, v ) such that v G L(W) 
and u =g v. But then, if Ui,u 2 G A* are geodesies with /(u! -1 ^) < 1 and 
ui,U2 reduce to v\,v 2 G L(W), then (ui, Vi), (u?, V2) G L(WD n ) whereas 
(vx,V2) G L(Dm)- This would show that the word-differences arising from 
{u, v) have bounded length, and so G is strongly geodesically automatic by 
Theorem 2.3.5 of [[Epstein et al., 1992|| . Thus it suffices to establish the claim 



that L(GW n ) is equal to the set of all geodesic words. 

Suppose that the claim is false, and let u be a minimal length geodesic 
word not contained in L{GW n ). Since L(GW n ) contains the empty word, 
u = wx for some word w and some x G A. Then w is a geodesic word and 
l(w) = l{u) — 1. Suppose that w reduces to v G L{W). Then w G L{GW n ) 
and (w,v) G L{GE n ). So (w,v) G L(WD n ). Therefore (u, v x) = (wx,vx) G 
L(WD n ). 

If vx G L(W), then by definition of GE n we get (u, vx) G L(GE n ), and 
so u G L(GW n ), a contradiction. On the other hand, if vx ^ L(W) and vx 
reduces to v' G L(W), then (v , v ') G L(M X ). From the definition of WDi we 
have (v x, v') G L{GE{), which implies vx G L(GW r ) for all r > 1. But then 
u G L(T n ), contrary to assumption. □ 

The above procedure is very similar to that described in [|Wakefield, 1997 . 



The principal difference is that our definition of the test-machines T n is rather 
simpler. Furthermore, in our implementation, we do not construct the non- 
deterministic automaton resulting from the quantification over the second 
variable in the definition of L(T n ). Instead, we construct the two- variable 
automaton with language 

{(w,u) I (w,u) G L(WD n ), u g L(GW n ), l(u) = l(w)}, 

and check during the course of the construction whether there are any words 
w ^ L(GW n ) that arise. If the construction completes and we find no such 
w, then L(T n ) is empty. Otherwise we abort after finding some fixed number 
of words (such as 500) w ^ L{GW n ) and use all of these words with their 
short-lex reductions to generate new word-differences. 

So the procedure stops if and only if G is word-hyperbolic. If it stops 
at the n-th stage, then WT> n is a finite set of word- differences which gives 
the best possible value of Papasoglu's constant for the particular choice of 
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generators. More precisely, the procedures described above find the set of all 
word- differences for all pairs of geodesies (w, u) which start at the same vertex 
of the Cayley graph and end at vertices at distance at most one apart, and 
where the word u is short-lex minimal. But by using a standard composite 
operation on two-variable automata as defined in [[Holt, 1996|| for example, we 
can easily compute from this the word- differences for general geodesic bigons 
in which at least one of the vertices is a vertex of the Cayley graph. Moreover, 
we have constructed an automaton whose language is the set of all geodesies 
in the Cayley graph that start and end at vertices of the Cayley graph. 



4 Finding the constant of hyperbolicity 

Throughout this section, we assume that G = {X) is a word-hyperbolic group 
and that 5 > is a constant such that all geodesic triangles in T X (G) are 
5-thin. The aim is to devise a practical algorithm to find such a 5, which 
should of course be as small as possible. As before, we assume that we have 
already calculated the short-lex automatic structure for G with respect to X. 

4.1 The reverse of a finite state automaton 

Our procedure makes use of reversed automata, and so we start with a brief 
discussion of this topic. 

If w G A* is a word, then we denote the reversed word by w R . Let 
M be a finite state automaton with alphabet A. We want to form the re- 
versed automaton M R with language {w R \w G L(M)}. We can define a 
non-deterministic version NM R of M R , simply by reversing the arrows of all 
transitions, and interchanging the sets of initial and accepting states. Then 
we can build a deterministic machine M R with the same language in the 
standard way, by replacing the set £ = S(NM R ) of states of NM R with its 
power set V(T,), and, for T, T G V(S), defining a transition in M R with a 
label x from T to T, if T is the set of all states v of M from which there 
exists an arrow labelled x in M with target some r G T. If we think of x as 
a partial map p x from S(M) to itself, then the existence of an arrow in M R 
from T to T is equivalent to saying that T is the full inverse image of T in 
iS(M) under p x . 

The initial state of M R is the set of all accepting states of M and a state 
of M R is accepting whenever it contains an initial state of M. 

In practice, we do not need to use the complete power set We 
start with the set of accepting states of M as initial state of M R , and then 
construct the accessible states and transitions of M R as the orbit of the initial 
state under the action of A* . 
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Let 0"o(M) and <j$(M r ) be the initial states of M and M R respectively. 
Suppose that v,w G A*, a (M) v = r, and a (M R ) w = T C V(E). Then, 
from the construction above, we see that r G T if and only if vw R G L(M). 
We shall need this property below, and when we compute the reverse of an 
automaton we need to remember the subsets of £ that define the states of 
M R . (This means that we cannot minimise M R , but in practice this does 
not appear to be a problem because, at least for word-acceptors of automatic 
group, M R does not seem to have many more states than M.) 

4.2 Reduction to short-lex geodesic triangles 

The reader should now recall the notation for 5-thin hyperbolic triangles with 
vertices a,b,c and edges u,v,w in the Cayley graph T = T X (G) defined in 
Section |^. We shall call a geodesic triangle in T short-lex geodesic, if its 
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Figure 3. A triangle with short-lex sides, annotated as in the text. 

vertices are vertices of T and if the words A* corresponding to the edges of 
the triangle (which we shall also denote by u,v,w) all lie in L(W); that is, 
they are all short-lex minimal words. 

It is important to work with short-lex triangles, because in general there 
are far more geodesic triangles and consideration of all of these is likely to 
make an already difficult computational problem impossible. 

Our algorithms are designed to compute the minimal 5 G N such that 
all short-lex geodesic triangles are <5-thin. In fact they do considerably more 
than this, because they compute the set of word-differences that arise from 
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the word-pairs (u, v), where u and v are the words from the two sides of such 
a triangle that go from a vertex of the triangle as far as the meeting points 
on the two sides. 

From the following proposition, we derive a bound on the value of the 
thinness constant of hyperbolicity for arbitrary geodesic triangles. As before, 
let T>m be the set of word-differences associated with the multiplier automata 
M x in the short-lex automatic structure of G, and let 7 be the length of the 
longest element in T>m- Let 7' be the length of the longest element in the 
final stable set WV n of word-differences defined in Section |3|. Although the 
bound in the next proposition will in general be too large, the fact that 7 
and 7' are usually smaller than 5 means that it is likely only to be wrong by 
a small constant factor. 

Proposition 4.1 Suppose that all short-lex geodesic triangles in F are 5-thin. 
Then all geodesic triangles in T are (5 + 2(7 + 7') + 3) -thin. 

Proof: Let A be any geodesic triangle in P We fix a direction around the 
triangle which we call clockwise. We define a new triangle A', which is equal 
to A except that its vertices are (if necessary) moved clockwise around A to 
the nearest vertex of P Phus each vertex is moved through a distance less 
than one 1, and it follows easily that the meeting points of the triangle move 
through a distance less than 2. It is not difficult to see that if p and q are 
A-companions (see Section |2] for definition), then there exist A'-companions 
p' and q' with d(p,p') + d(q, q') < 2. (A careful argument is needed if p and 
q are near meeting points.) It follows that, if A' is <5'-thin for some 5', then 
A is (5' + 2)-thin. 

Let a, b, c be the vertices of A' in clockwise order. Let a' be the vertex on 
the union of the three sides of A', adjacent to a in an anticlockwise direction, 
and similarly for b' and c'. (So the vertices of A lie on the edges (a', a), (&', b) 
and (c', c).) We will assume that the six vertices {a, a', b, b', c, c} are distinct, 
and leave to the reader the minor modifications needed for the cases where 
two or more of them coincide. 

Let u, v and w be the paths from b to c, c to a and a to b, respectively, 
defining short-lex reduced words u, v and w in L(W), and let A" be the 
resulting short-lex geodesic triangle with vertices a, b and c. Phe triangle A' 
is not necessarily geodesic, but the paths u', v', w' on A' from b to c', c to a' 
and a to b', respectively, are geodesic, because they are subpaths of the sides 
of A. 

Let u", v" and w" be the paths from b to c', c to a' and a to b', respectively, 
defining short-lex reduced words in L(W). Phe situation is illustrated in 
Figure f|. 

Phen u and u" are words in L(W) that have the same starting point in F 
and end a distance one apart. So the word-differences arising from (u,u") lie 
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■ denotes a vertex of A 

Figure 4. This diagram illustrate the situation described in the proof of Propo- 
sition O. 
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in Dm- Since the paths u' and u" have the same starting and end points, vl 
is a geodesic and u" G L(W), the word-differences arising from (u',u") lie in 
WT> n . It follows that the word-differences arising from (u,u') have length at 
most 7 + 7', and similarly for (v,v') and (w, w'). 




Figure 5. These illustrates the origin of 5 + 2(7 + 7') + 1 in the expression 
bounding d(p',q'). In the diagram on the left, l(ca) = l(ca') + 1, and in the 
diagram on the right, /(ca) = l(ca'). The shapes of the curves are due to the 
fact that we are sometimes looking at points which are equidistant from a and 
sometimes at points which are equidistant from c. Interested readers are left to 
work out the details for themselves. 

The triangles A' and A" have the same vertices a, b and c. Each side of 
A' is at least as long as the corresponding side of A", but no more than one 
unit longer — this can be deduced from the fact that the sides of the original 
triangle A are geodesic. It follows that the distance from a vertex of A' to 
its two adjacent A'-meeting points is within one unit of its distance to its 
two adjacent A"-meeting points. We then see that, if all pairs p, q of A"- 
companions satisfy d(p, q) < 5, then all pairs p', q' of A'-companions satisfy 
d(p', q') < 5 + 2(7 + 7') + 1. The result now follows. 

Part of the argument is illustrated in Figure |5| □ 

4.3 The automata FRD and FRD 3 

In this section, we describe a finite state automaton which we shall call FRD, 
which stands for 'forward, reverse, difference'. Roughly speaking, it is a two- 
variable machine which reads the two sides emerging from a vertex of a short- 
lex geodesic triangle as far as the meeting points on those two sides. We also 
describe an associated automaton FRD 3 which consists of three copies of 
FRD. The three pairs of words read by FRD 3 will be accepted when they 
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are the three pairs of edges emerging from the three vertices of some short-lex 
geodesic triangle, ending and meeting at the meeting points of the triangle. 




Figure 6. This diagram shows the meeting vertices, denoted by d, e and /. 
Companions have been joined by a line. The diagram on the left shows the 
situation where the perimeter of the triangle is even and the diagram on the right 
where the perimeter is odd. 

When the perimeter l(u) + l(v) + l(w) of a short-lex geodesic triangle is 
even, the meeting points d, e, / are vertices of T which lie on u = be, v = ca 
and w = ab respectively. When the perimeter is odd, however, each meeting 
point lies in the middle of an edge of T, which is not convenient for us. We 
therefore move them to a neighbouring vertex and re-define d, e, / to be the 
points on u,v,w at distance p(b) + 1/2, p(c) + 1/2 and p(a) + 1/2 from b, c 
and a, respectively, and call d, e, / the meeting vertices. Each of the three 
vertices of the triangle has two incident sides, one in the clockwise direction 
and the other in the anti-clockwise direction, starting from the vertex. In 
these terms, if we start at a vertex of the triangle and move away from this 
vertex along the two sides of the triangle emerging from it, one edge at a 
time, then we have to move one extra edge along the clockwise side than the 
anti-clockwise side to reach the meeting vertices on the two sides. 

Let W be the word acceptor for the short-lex automatic structure of G. 
We assume that the reverse W R of W has been computed, as described in f4TT| . 
For a short-lex geodesic triangle with meeting vertices d, e and / defined as 
above, we denote the elements of G corresponding to paths from d to /, e to 
d, and / to e, by r, s and t, respectively. This is illustrated in Figure |[ 

Let T>\ be the set of elements r of G that arise by considering all such 
triangles. (By symmetry, this set includes all of the elements s and t as well.) 
Let T>2 denote the set of all elements of G of the form 

{w(i)~V(2) | i e Z, < % < p{a)}, 

for all triangles under consideration; that is, the set of word-differences arising 
from reading the two edges of the triangle from a up to the meeting points. 
(For triangles with even perimeter, the elements r, s, t lie in both T>\ and T>2-) 
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Let T>t — Di U T>2- Then our assumption that geodesic triangles are <5-thin 
implies that T>t is finite, and that its elements have length at most 5+1. 

The automaton FRD is defined as follows. Its states are triples (a, E, g>) 
with a G S(W), E G S(W R ) and g G V T . (This explains the name FRD.) 
For x G A and y G A\ a transition is defined with label (x,y) from (a, E, g) 
to (a', E', (/) if and only if o"^ = a', TP = E' and x~ x gy =q g' ■ (We allow y 
but not a; to be the padding symbol because, in the case of triangles with odd 
perimeter, we need to read one extra generator from the left hand edge than 
the right hand edge to get to the meeting vertex. We define E $ = E for all 
S.) The initial state is (o- Q (W), a (W R ), 1 G ). The accepting states of FRD 
are those with the third component g G T)\. 

The above description is not quite correct. In fact, each state has a fourth 
component which is either or 1, and is 1 only when a pair (x, $) has been 
read. There are no transitions from such states. 

We also need to consider an automaton FRD 3 , which consists of the prod- 
uct of three independent copies FRD a , FRDb and FRD C of FRD. Its input 
consists of sextuples of words (u a ,v a ,Ub,Vb,u c ,v c ) G (A*) 6 , where iu a ,v a ), 
(ub,Vb) and (u c ,v c ) are input to the three copies of FRD. A state of FRD 3 
consists of a triple (r a , u, r c ), where r a = (a a , E a , g a ) is a state of FRD a , etc. 
The initial state of FRD 3 is the triple consisting of the initial states of the 
three copies. 

To specify the set A(FRD 3 ) of accepting states, we recall that a state 
E of W R is a subset of the set S(W) of states of W. We have (r a ,r b ,r c ) G 
A(FRD 3 ) if and only if a a G S 6 , cr b G S c , a c G E , and g c g b g a = 1 G . 

The proof of the following lemma should now be clear. 

Lemma 4.2 The triple (u a , v a , Ub, Vb, u c , v c ) is accepted by FRD 3 if and only 
ifubV R , u c v R and u a v R all lie in L(W) and form the three sides of a short-lex 
geodesic triangle in T. 

4.4 Proving correctness of FRD 

If we can construct FRD, then we can compute the value of the hyperbolic 
thinness constant 5 for short-lex geodesic triangles as the maximum of the 
lengths of the words in T>t- (Or more exactly, this number plus 1, because of 
our re-definition of the meeting points of the triangles.) Conversely, we have 
already computed W and W R , so the construction of FRD only requires us 
to find the set V T . The idea is to construct a candidate for FRD, and then 
to verify that it is correct. 

In one of our short-lex geodesic triangles, two of the edge- words u, v can 
be chosen as arbitrary elements of L(W), and the third is then determined 
as the representative in L(W) of (uv) -1 . So we can proceed by choosing 
a large number of random pairs of words in u,v G L{W) (for example we 
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might choose 10000 such pairs of length up to 50), computing w as described, 
and then computing the set T>t of word-differences arising from all of these 
triangles. We do this several times until the set V T appears to have stopped 
growing in size. (We are assuming that G is word-hyperbolic, so we know 
that the correct set V>t is finite.) We then proceed to the verification stage, 
which is computationally the most intensive. 

The idea is to compute a two- variable finite state automaton GP (geodesic 
pairs) of which the language is the subset of A* x A* defined by the expression 

{( Wl ,w 2 y EA* x A*\ 

3(tt , v a , u b , v b , u c , v c ) e L(FRD 3 ) : w x = u a v R , w 2 = v a uf}. 

Then GP accepts the set of pairs of sides (w, v R ) emerging from the vertex a 
in the triangles that are accepted by the current version of FRD 3 . (During 
the course of the program, FRD changes as more information is incorporated 
into it.) Thus FRD is correct if and only if L(GP) = L(W) x L(W R ). Since 
checking for equality of the languages of minimised deterministic automata 
is easy, we can perform this check provided that we can construct GP. 

Furthermore, if the check fails then our definition of GP ensures that 
L(GP) C L(W) x L(W R ). So we can find one or more specific words 
(wi,w 2 ) £ L(W) x L(W R ) \ L{GP) and then compute the word-differences 
arising from the short-lex geodesic triangle having wi and w R as two of its 
sides. We can then adjoin these to T>t and return to the construction of 
FRD. 

The construction of GP can be carried out in principle, but because of 
the large number of quantified variables involved in the above expression, a 
naive implementation would be hopelessly expensive in memory usage. 

We now give a second and more detailed version of our implementation 
of this construction, in a way that makes the computation easier to carry 
out. It remains heavy in its memory usage, but it does at least work for easy 
examples. The basic objects collected during the course of the computation 
are the word- differences T>t referred to above, which are used in constructing 
FRD, and the triples (r, s,t) of small triangles in the Cayley graph whose 
vertices are the meeting vertices of some short-lex geodesic triangles. The 
situation is shown in Figure |]. 

The general idea is to map out a short-lex geodesic triangle by advancing 
from one vertex, say a, using FRD. Then, at a certain moment, the machine 
'does the splits', with a non-deterministic jump corresponding to one of the 
small (r, s, t)-triangles. To complete the triangle, the two legs have to be 
followed to the other vertices b and c, respectively, of the large triangle. Each 
of the legs follows the reverse of FRD. 

Here is a third, even more explicit, version of the construction, which can 
be skipped by readers who are only interested in the conceptual description 
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of the program. We assume that we have constructed a candidate for FRD 
explicitly. We do not construct FRD 3 explicitly, but we do make a list of all 
triples (r,s,t) such that r,s,t G S(FRD) and (r,s,t) G A(FRD 3 ). Having 
done that, we can forget the structure of the states of FRD as triples, and 
simply manipulate them as integers. 

We also need a version of an automaton FRD R that accepts the reverse 
language of FRD. In this case it is convenient to work with a partially non- 
deterministic version — that is, it is deterministic except that there are many 
initial states, not just one. The states are subsets of S(FRD) as described 



in |4.1| and the transitions and accepting states are also as described there. 
But instead of having a unique initial state, for each accepting state r of 
FRD we make the singleton subset {r} into an initial state of FRD R . Note 
also that if (u,v)^ G L(FRD) with l(u) = l(v) + 1, then the reversed pair 
accepted by FRD R is (u R , $v R ); that is, the padding symbol comes at the 
beginning rather than the end. In general, let us denote the word-pair formed 
by inserting the padding symbol at the beginning by '(u,v). 

We shall now describe a non-deterministic version NGP of GP (geodesic 
pairs). Subsequent to its construction, it can be determinised using the usual 
subset construction, minimised, and its language compared with L(W) x 
L{W) R . 

The precise description of NGP is rather technical, so we shall first at- 
tempt to explain its operation. A (padded) pair of words (wi,w 2 y is to be 
accepted if and only it satisfies the expression above; that is, (wi,W2) = 
R ,v a u R ), where (u a ,v a y G L(FRD), and there exist words u b ,v c G A 
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with (u b , v b ) ] , (u c , v c y G L(FRD) and (u a ,v a ,u b ,v b ,u c ,v c ) G L(FRD 3 ) 

Equivalently, (writing the reverses of u b , u c , v b and v c by capitalising u 
and v), (wijW^y is accepted if and only of (w\,W2) = (u a V b ,v a U c ) where 
(u a ,v a y G L(FRD), and there exist words U b , V c G A* with ^(U b ,V b ), 
] (U C ,V C ) G L(FRD R ) and 



(u a , v a , U R , V b , U R , V R ) G L(FRD 3 ). 

So the accepting path of (11)1,11)2) through NGP will be in two parts, the 
first (u a ,v a ) and the second (V b , U c ). A picture of this is shown in Figure [7|. 
Furthermore, we have either l(u a ) = l(v a ) or l(u a ) = l(v a ) + 1, depending on 
whether the perimeter of the geodesic triangle which has W\ and W2 as two of 
its sides is even or odd. In the first case, it is possible for (V b , U c ) to be empty, 
which occurs when the vertices b and c of the geodesic triangle coincide. 

The first part of the path through NGP is simply a path through FRD, 
ending at a state a G A(FRD). In Figure 0, an intermediate state is also 
denoted by a. 

The second part corresponds to the two paths (U b , V b ) and (Ub, U c ) through 
FRD R . These paths must end at the unique accepting state of FRD R . This 
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part of NGP is non-deterministic, because we need to quantify over their 
second variables as described in Section 0. The initial states it\ = {ai} and 
7r 2 = {cr 2 } of these two paths through FRD R must be such that (a, ai, <r 2 ) G 
A(FRD 3 ). This is equivalent to (u a , u OJ U?, V b R , U R , V R ) G L(FRD 3 ). In 
Figure [7| intermediate states are denoted by (px,p 2 ). 

In our implementation of NGP, we prefer to avoid e-transitions, and so 
the non-deterministic jump from the first to the second part of the path is 
combined with the first transitions in the second part of the path. In the case 
where there is a padding symbol, the last transition in the first part of the 
path is combined with the first transitions in the second part. An advantage 
of this is that we can eliminate the use of the padding symbol in the middle 
of a word, which can otherwise be quite troublesome to deal with (in terms of 
writing special code to take the unnecessary padding symbols into account). 

The jump also introduces a large amount of non-determinism into NGP. 

The states of NGP are triples (a, pi,p 2 ), where a G S(FRD) U {00} and 
p±, pi G S(FRD R ) U {0,oo}. For each such state either p 1 = p 2 = and 
a ^ 00. or a = 00 and p 1 7^ 7^ p 2 - Informally, and 00 as just introduced 
have the following significance. In the course of accepting a string, the three 
components a, p\ and p 2 each have to pass through FRD exactly once. More 
precisely, the component o passes once through FRD during the first part of 
the path and the components p\ and p 2 pass through FRD R once during the 




First part shaded 



a 

Figure 7. This shows the path in the automaton, broken into a first part which is 
in the automaton FRD and a second part which is in two copies of the automaton 

FRD R . 
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second part of the path, each time moving from an initial state to an accept 
state of FRD or FRD R respectively. It is convenient to assign the name oo 
to the state a after it has completed its passage through FRD. We need to 
attach names to the states pi (i = 1, 2) during the first part of the path, and 
we attach the name to remind us that p$ has not yet started its passage 
through FRD R . When a padding symbol is read in w iy the state pi is set to 
oo to remind us that the passage of pi through FRD R is now complete. In 
other words, the state is set to oo the next move after arriving at the vertex 
of the triangle. We do not allow p\ = p 2 = oo, because we will stop if we 
reach both vertices b and c simultaneously. 

We can save space by storing the triple (a, 0,0) as a pair (cr, 0) and the 
triple (oo,pi,p 2 ) as a pair (pi,p 2 ). It is easy to see that this captures all the 
information. In this discussion, we continue to use more revealing triples, 
rather than more concise pairs. 

The unique initial state is (a (FRD), 0, 0). The accepting states are 
(oo,pi,p 2 ) where either p 1 = oo or pi lies in A(FRD R ) and the same is 
true for p 2 . The reader may like to be reminded that a state of A(FRD R ) 
is an accept state of FRD R , and that this is a subset of the set of states of 
FRD which contains the initial state ao(FRD). 

There is also another kind of accept state, corresponding to the situa- 
tion b = c, that is, that there are two different geodesies from a to b = c. 
This means that W\ =g w 2 . Since we are dealing with short-lex geodesies, 
Wi will be short-lex from a to b and w 2 will be the reverse of a short- 
lex geodesic from c to a. Such an accept state has the form (cr, 0) where 
(cr, a (FRD), a (FRD)) G A(FRD 3 ). 

There are other degenerate situations, but the others are all covered by 
the main description, as the reader can easily verify. 

There are three types of transitions of NGP, which we shall now de- 
scribe. In general, we denote the label of such a transition by (x a b,x ac ), 
where x a b,x ac G (The idea is that x a b represents a variable generator 
in the path from a to b.) The first two types of transition correspond to 
the transitions of the first and second parts of the accepting path of (wi, w 2 ) 
through NGP, and the third type of transition corresponds to a jump from 
the first to the second part. 

Transitions of the first type have the form (cr, 0, 0) — > (r, 0, 0) where there 
is a transition o — > r of FRD with the same label. However, we must have 
x a b, x ac G A; that is, x ac is not allowed to be the padding symbol $. (Any 
such transition (x, $) of FRD will be absorbed into the jump, and combined 
with the first transitions of the two copies of FRD R after the jump, which 
are bound to have labels of the form (y, $) as we see from Figure [| Recall 
that we do not want the padding symbol to occur in the middle of either of 
the words w±, w 2 . The strategy explained here avoids that danger.) 

Transitions of the second type have the form (oo,7ri,7r 2 ) — > (oc,pi,p 2 ). 
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They occur whenever there exist Xdb, Xdc G and transitions tt% — > p\ and 
7r 2 — > P2 of FRD R with labels (a;^, x a h) and (x ac ,Xrf c ), respectively. There is 
the further restriction that x a b = $ if and only if Xdb = $• This means that 
our path has previously arrived at the vertex b — see Figure [5]. In this case we 
have pi = oo and either tti = oo or n 1 G A(FRD R ). Similarly, x ac = $ if and 
only if Xdc = $• In this case p 2 = oo and either tt 2 = oo or n 2 G A(FRD R ). 
As explained above, we cannot have x a b = x ac = $. In other words, padding 
symbols occur only at the end of at most one of the words Wi,w 2 . 

The transitions of the third type are jumps from (a, 0, 0) to (oo,p 1; p 2 )- 
These are of two subtypes, depending on whether the geodesic triangle defined 
by the accepting paths that pass through them has perimeter of even or odd 
length. 

Those of even perimeter subtype occur when there exist Xdb, %dc £ A 
and initial states 7i"i = {o~i}, 7t 2 = {02} of FRD R with the property that 
(a, <Ji,a 2 ) G A(FRD 3 ). Furthermore there are transitions ni — > pi and 
^2 - ► P2 with labels (x^,x a b) and (xac^dc) respectively. There is the further 
restriction that x a b = $ if and only if x,j,b = $, and in this case we have p\ = 00 
and 7Ti G A(FRD R ). Similarly, x ac = $ if and only if Xdc = $, and in this 
case P2 = 00 and 7r 2 G A(FRD r ). 

Those of the odd perimeter subtype arise only for x a b,x ac £ A, and they 
occur when there is a transition a — > cr' with label (x a &, $) of FRD. Further- 
more, there exists 2;^ G A and initial states 7i"i = {o"i}, tt2 = {o"2} of FRD R 
with the property that (a', <Ti, 02) G A(FRD 3 ). Also, there are transitions 
7Ti — > pi and 7r 2 — > P2 with labels (xdb, $) and (x ac , $) respectively. 

5 Examples 

In this final section, we describe the performance of these algorithms on the 
following four examples. 

Gi = (a,b,c,d\a~ 1 b~ 1 abc~ 1 d~ 1 cd = l), 
G 2 = (a,b\a 2 = b 3 = (ab) 7 = 1), 
G 3 = (a,b\ {b- x a 2 ba- 3 f = 1) and 

G 4 = (a, b, c, d, e, / | a 4 = b 4 = c 4 = d 4 = e 4 = f 4 = 

aba~ x e = bcb~ 4 f = cdc~ 4 a = ded~ l b = efe~c = faf~ 4 d = 1). 

Of these, G% is a surface group of a genus two torus, G 2 is the (2, 3, 7)- 
von Dyck group, G3 is obtained from one of the well-known family of non- 
Hopfian Baumslag-Solitar groups by squaring the single relator, and G 4 is 
the symmetry group of a certain tessellation by dodecahedra of hyperbolic 
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3-space as featured in the video 'Not Knot' ( ||Not Knotf ). They are all word- 
hyperbolic groups. 

For all four the verification of hyperbolicity, as described in Section |3], was 
relatively easy, with the first three examples completing in a few seconds, 
and G 4 requiring about 20 seconds cpu-time. We present details of these 
calculations in Table 1. The automata WD n , GE n , GW n and T n 
described in Section HL and the constants 7 and 7' are as defined before 



Proposition |4.1| . The notation m — > n in the table means that an automaton 
had m states when it was first constructed, and it was then minimised to 
an equivalent automaton with n states. The last example demonstrates the 
phenomenon that the automata involved are smaller when the data is correct. 



Table 1. Verifying Hyperbolicity 



Grp 


n 


S(WD n ) 


S{GE n ) 


S(GW n ) 


S(T n ) 


7 


7' 


G l 


1 


33 


121 -> 49 


49 -> 49 


265 


4 


4 


G 2 


1 


30 


627 -> 92 


80 -> 52 


936 








2 


32 


664 -> 94 


78 -> 54 


769 


7 


7 


G 3 


1 


55 


689 -> 136 


152 -> 96 


1270 


6 


6 


G4 


1 


75 


896 -> 284 


454 -> 409 


10635 








2 


97 


1135 -> 309 


443 -> 378 


12407 








3 


103 


1211 -> 318 


424 -> 63 


1713 


4 


4 



In Table 2, we present details of the calculation of the thinness constant 
for short-lex geodesic hyperbolic triangles in the first three of the examples. 
The set V and the automata FRD, FRD 3 , NGP and GP are as defined in 
Section ^ (where NGP is the non-deterministic version of GP). The separate 
lines of data for each group represent successive attempts at the computation, 
with the last line representing the correct data. After each attempt, the 
automaton with language L = L(W) x L(W R ) \ L{GP) was constructed 
and, when it was nonempty, words (w\, W2) G L were found and used to find 
an improved set T>. The language L was found to be empty after the final 
computation for each group, thereby proving correctness of the data. 

The behaviour of G3, which is the most difficult example for which we 
have successfully completed the calculations, is probably the best indicator 
of the way in which more difficult examples are likely to behave. For exam- 
ple, the largest and most memory intensive part of the computation is the 
determinisation of NGP, and many parts of the calculations are significantly 
more expensive on the earlier passes, when the data is incorrect, than in the 
final correct stage. These computations were carried out using a maximum 
of 256 megabytes of core memory and about the same amount of swap space. 

We have not yet been able to complete the calculations for G4; indeed we 
have not progressed further than the first construction of NGP, which has 
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several million states. We will need more memory (probably more than a 
gigabyte) if we are to proceed to construct the determinised version GP. 



Table 2. Finding the Constant of Hyperbolicity 



Grp 


V 


S(FRD) 


A(FRD' A ) 


S(NGP) 


S(GP) 


5 


G 1 


25 


137 


65785 


12249 


8049 -> 2185 






49 


169 


65857 


12281 


5457 -> 625 


4 


G 2 


104 


1174 


73822 


89802 


35824 -> 4904 






111 


1199 


74047 


90450 


31374 -> 1508 


7 


G 3 


71 


755 


795436 


274186 


1872679 -> 531434 






195 


1430 


801745 


280240 


1695944 -> 443570 






241 


1741 


806923 


284328 


1237158 -> 85044 






257 


1845 


807136 


284478 


676645 -> 3803 
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